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Abstract 

Computational Fluid Dynamics based de- 
sign methods are maturing to the point that 
they are beginning to be used in the aircraft 
design process. Many design methods however 
have demonstrated deficiencies in the leading 
edge region of airfoil sections. The objective 
of the present research is to develop an efficient 
inverse design method which is valid in the lead- 
ing edge region. The new design method is a 
streamline curvature method, and a new tech- 
nique is presented for modeling the variation of 
the streamline curvature normal to the surface. 
The new design method allows the surface co- 
ordinates to move normal to the surface, and 
has been incorporated into the Constrained Di- 
rect Iterative Surface Curvature (CDISC) de- 
sign method. The accuracy and efficiency of 
the design method is demonstrated using both 
two-dimensional and three-dimensional design 
cases. 

Nomenclature 

b wing span 

c airfoil or wing local chord 

C curvature 

C p surface static pressure coefficient 

M Mach number 

Re Reynolds number based on chord 

s arclength 

V velocity tangent to streamline 
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x, y, z Cartesian coordinates 
y + incompressible law-of-the-wall coordinate 

a angle-of-attack, degrees 

r] direction normal to surface 

a scaling coefficient 

Subscripts 

L lower surface 

U upper surface 

o surface value 

oo free stream value 

Superscripts 


* sonic value 

Introduction 

Computational Fluid Dynamics (CFD) 
methods have matured to the point that they 
are being used to analyze full aircraft configu- 
rations [1], as well as to develop experimental 
testing techniques [2]. One area of particular 
interest to the aircraft industry is the use of 
state-of-the-art CFD methods in the design pro- 
cess. This interest has led to the development of 
many efficient design methods which can be di- 
rectly coupled with Navier-Stokes solvers. The 
CFD based design methods must be efficient 
and accurately model relevant flow physics. One 
such example is the Constrained Direct Iterative 
Surface Curvature (CDISC) method developed 
at the NASA Langley Research Center [3]. 
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The C'DISC design method uses analyti- 
cally derived pressure/geometry relationships to 
modify a given surface geometry to match a 
prescribed target pressure distribution. The 
method is valid for subsonic, transonic, and su- 
personic flow regimes. The target pressure dis- 
tribution can be specified by a designer based 
on their own interests. For example, the target 
pressure distribution may be specified to give 
a favorable pressure gradient over a portion of 
an airfoil [4,5]. The C'DISC method also offers 
many geometric constraints which can be im- 
posed during a design study. A complete de- 
scription of the method can be found in Refer- 
ence 3. 

Design methods such as C'DISC have tra- 
ditionally been used to improve the cruise per- 
formance of a configuration. As these design 
methods have matured, interest has grown in 
applying such design methods to improve the 
configuration in other important flight regimes. 
One area of significant interest is the design and 
optimization of high-lift systems deployed dur- 
ing takeoff and landing [6]. For such design 
cases, surface geometry modifications are gen- 
erally limited to small regions, such as the lead- 
ing edge region of a flap element. Thus a design 
method must work well in leading edge regions 
before it can be applied to such challenging de- 
sign tasks. 

Many design methods including C'DISC' 
have demonstrated some deficiencies in the lead- 
ing edge region [5,6]. This is not surprising 
given the strong flow gradients and high surface 
curvature in this region. Several approaches 
have been developed to improve the leading 
edge region, which include representing the sur- 
face geometry using families of smooth analytic 
functions [6,7]. Such techniques can be useful, 
but their use can raise concern of biasing a de- 
sign by restricting it to fit within a certain geo- 
metric family. 

The objective of the present research is to 
develop a design method which performs well in 
the leading edge region at subsonic speeds. The 
method is based upon the same streamline cur- 


vature principles as the C'DISC' design method, 
however a new method is presented for mod- 
eling the variation of the streamline curvature 
normal to the surface. The new design method 
allows the surface coordinates to move normal 
to the surface. The new design method, SC'2D, 
has been incorporated into the C'DISC' design 
method. Several two-dimensional and three- 
dimensional design cases are presented which 
demonstrate the efficiency and accuracy of the 
SC'2D design method. The new design method 
will be shown to perform quite well in the lead- 
ing edge region. 



Figure 1. Flowchart for the C'DISC' design 
method. 

CDISC Design Method 

A flowchart for the CDISC design method is 
shown in Figure 1. The CDISC module is indi- 
rectly coupled to a flow solver using a script hie. 
After a prescribed number of how solver itera- 
tions the current grid and solution restart hies 
are passed to the CDISC module. The surface 
geometry and pressure distributions are first ex- 
tracted from the grid and restart hies. The 
target pressures are modified to meet any pre- 
scribed how constraints, then passed to the ba- 
sic surface design module. This is where the new 
SC'2D design method has been added, which 
will be discussed in detail below. After the 
basic surface design is complete, the geometry 
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is modified to meet any geometry constraints. 
The volume grid is then updated and returned 
to the flow solver for analysis. This process is 
defined as a design cycle. This procedure con- 
tinues through several design cycles until the 
design is within acceptable agreement with the 
specified targets. 

SC2D Design Method 

The new design algorithm is based upon 
the same streamline curvature relationship used 
by Barger and Brooks [8], and more recently 
used by Campbell [3]. The relationship is de- 
rived from the momentum equations written in 
a streamline coordinate system; 


curvature near the surface. Figure 2 shows a 
streamline pattern near the surface of a two- 
dimensional body. Moving away from the sur- 
face a small distance, the streamline curvature 
is nearly equal to the surface curvature. The 
streamline curvature is thus assumed to vary in 



— = —C(i])(h] (f) 

where V is the velocity tangent to the stream- 
line, C is the streamline curvature, and i] is 
the direction normal to the streamline. Before 
equation (f) can be used to develop a design 
algorithm, the variation of curvature normal to 
the streamline must be formulated. Barger and 
Brooks [8] assumed that the streamline curva- 
ture decays exponentially into the free stream; 

C{rj) = C 0 e ~ k " (2) 

where C 0 is the curvature at the airfoil surface, 
and k is a constant. Downstream of an airfoil 
leading edge this appears to be a reasonable as- 
sumption, and has been used successfully for 
designing transonic airfoils and wings [4]. In 
the leading edge region however, Barger and 
Brooks noted that this expression may not be 
valid. In this region the streamline curvature 
along a normal to the surface may change sign, 
particularly near the stagnation point. Equa- 
tion (2) does not account for a change in the 
sign of the streamline curvature. As a result, a 
design method based on equation (2) may be- 
come problematic in the leading edge region [5]. 

In the present research a new approach 
is presented for modeling the variation of the 
streamline curvature normal to the surface. For 
an iterative design procedure, it is only neces- 
sary to consider the variation of the streamline 


Figure 2. Near surface streamline pattern. 

a linear manner near the surface; 

C{i]) = C 0 + di] (3) 

where d is a constant. This form is valid in 
the leading edge region, and does allow for a 
change in the sign of the streamline curvature. 
This expression is substituted into equation (I), 
and integrated from the surface (?/ = 0) to the 
point Pi near the surface (?/ = ?q). This gives; 

1,1 (fe) = “ + 1"' 2 ) 141 

where V 0 is the inviscid surface velocity. Since 
the point Pi is a small distance from the surface, 
small perturbations can be used; 

Id = V 0 + AV 
m = do + a i] (5) 

which are used to give: 

In (l + = - (c 0 Ari + . (6) 

Consistent with the small perturbation ap- 
proach, the natural logarithm can be written 
as a Taylor series. Discarding the higher order 
terms yields second order accuracy: 

« -Co A t/ Vo / 0. (7) 


3 



Equation (7) relates the difference in the 
velocity between two streamlines to the nor- 
mal distance between them. The inviscid ve- 
locity in equation (7) can be replaced with the 
inviscid Mach number by assuming a zero tem- 
perature gradient normal to the surface. The 
inviscid Mach number is computed from the 
surface pressure coefficient using isentropic re- 
lations. The surface normal perturbation can 
now be solved for: 


A?/ 0 ~ 


AM 1 


( 8 ) 



Figure 3. Leading edge region of an airfoil 
section. 


Equation (8) has two singularities; M 0 = 0* 
and Co = 0. The stagnation point is where M 0 
= 0, and is avoided implicitly during the design 
procedure as discussed below. The denomina- 
tor is modified to avoid zero curvature. Lastly 
a scaling coefficient is added to control the de- 
sign procedure and maintain numerical stabil- 
ity. The surface normal perturbation relation is 
thus; 


A?/ 0 



1 

1+ | C 0 


( 9 ) 


where a is the scaling coefficient, and is gener- 
ally less than 0.40. 


Surface Design Procedure 


The steps used to design an airfoil section 
using equation (9) are now presented. Figure 
3 shows a sketch of an airfoil section leading 
edge with both the stagnation streamline and 
stagnation point shown. The stagnation point 
serves as the delimiter between the upper and 
lower surfaces of the airfoil. This is in contrast 
to many design methods which begin at the ge- 
ometric leading edge point. The airfoil is pa- 
rameterized by arclength, with the stagnation 
point as the origin. The curvature for a given 
surface is defined as: 


C(s) 


^S^SS Zs%ss 
( Xs 2 + ~~ s 2 ) 3/2 


( 10 ) 


where the subscripts denote differentiation with 
respect to arclength. This formulation avoids 
singularities in the leading edge which can occur 
if the curvature is parameterized by x. 


The next step in the design process is to 
compute the perturbation of the surface nor- 
mals along the airfoil. It is important to note 
that the stagnation point is treated as a fixed 
point during this step, avoiding the singularity 
of M 0 = 0 in equation (9). At each point along 
the airfoil the surface point is moved along the 
surface normal by the amount A?/ 0 computed 
using equation (9). This process is sketched 
in Figure 4 for point i along the surface. This 



s 



along local surface normal. 


method allows both the x and .r coordinates to 
move, with the upper and lower surfaces of the 
airfoil generally rotating about the f ixed stag- 
nation point. The leading edge point generally 
moves, and a search is performed to find the 
leading edge point of the new airfoil. It should 
also be noted that the upper and lower surface 
trailing edge points are allowed to move inde- 
pendently during the design. This typically al- 
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lows a blunt trailing edge to develop, which may 
not be desirable for all design cases. If a closed 
trailing edge is desired, the upper and lower sur- 
faces of the new airfoil are rotated about the 
new leading edge point to close the trailing edge. 

The ability of the design method to move 
both the x and .r coordinates allows the chord 
length to vary during a design if desired. For 
the current research however, the chord length 
is conserved. The new airfoil coordinates are 
splined versus the original airfoil arclength dis- 
tribution to maintain a uniform arclength dis- 
tribution during the design. The final step is 
the explicit smoothing of the new surface coor- 
dinates. Typically one smoothing pass is used 
per design cycle using the existing smoothing 
method in CDISC [3]. 


techniques are used to accelerate the conver- 
gence to steady state. Additional details of the 
code are found in Reference 9. The adequacy of 
the code for examining the low-speed flows in 
the current study is presented in Reference 10. 

Initial and Target Airfoil Sections 

The airfoil sections used to develop the vari- 
ous test cases for the present research are shown 
in Figure 5. The NAC'A 0012 airfoil section 
was used to formulate the initial converged flow 
held solutions for the various test cases. The 
remaining airfoil sections were used to gener- 
ate converged how held solutions which provide 
the “target” pressure distributions for the de- 
sign test cases, which will be discussed below. 


Transient Leading Edge Treatment 

Experience with the SC2D design method 
has demonstrated that the leading edge region 
can still be problematic during the early de- 
sign cycles. To address this “start-up” problem, 
a transient polynomial patching technique was 
developed. During the early design cycles, a 6th 
order polynomial is htted to the leading edge re- 
gion. This typically replaces only the first 0.5 - 
1.0% of the airfoil. As the lift coefficient nears 
convergence this patch is removed from the de- 
sign procedure, allowing the airfoil to develop 
in a non-constrained manner. 



Figure 5. Initial and target airfoil sections. 


Results 


Flow Solver 

The flow solver used for the Euler and 
Navier-Stokes simulations is TLNS3D-MB. 
The code solves the time-dependent, three- 
dimensional, thin-layer compressible Navier- 
Stokes equations on block-structured, body- 
fitted grids. The equations are discretized in a 
central difference finite volume formulation, and 
integrated using an explicit second-order accu- 
rate Runge-Ivutta time-stepping scheme. Multi 
grid, grid sequencing, and local time stepping 


The accuracy and efficiency of the new 
SC2D design method will be demonstrated with 
several two-dimensional and three-dimensional 
test cases. The Euler and Navier-Stokes design 
cases were performed on a Cray C-90. The com- 
putational times for the target solutions pre- 
sented below are based on the lift coefficient 
converging to four decimal places. The panel 
method results were run on an SGI Octane with 
an R12000 processor. For this single test case 
the SC2D design method was not run within the 
CDISC design method. 
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Inviscid Airfoil Design 

The first test case is an inviscid airfoil design 
for M (:>0 = 0.30 and a = 4.00°. The target air- 
foil is the RAE 5212 which was chosen due to its 
aft camber and non-parametric nose shape. Fig- 
ure 6 examines the influence of grid refinement 
on the RAE 5212 surface pressure distribution. 
The grid dimensions such as 289 X 33 repre- 
sent the number of grid points in the stream 
wise and normal directions, respectively, of the 
C'-grid topology. The results for both grids are 
nearly identical, and the fine grid was used for 
this test case. 



X/C 


Figure 6. Influence of grid refinement 
on the RAE 5212 airfoil, inviscid solution 
(Moo = 0.30, a = 4.00°). 

The results of this design test case are pre- 
sented in Figure 7 where the initial and target 
results are shown. The .r airfoil coordinates are 
plotted on an expanded scale for clarity. The 
SC2D design method has done an excellent job 
of replicating the target pressure distribution in 
only 25 design cycles. The maximum difference 
between the target and design pressure distribu- 
tions is less than ± 0.005 over the entire airfoil. 
The design converged rapidly in a monotonic 
fashion, with the lift coefficient converging to 
within 98% of the final value after only 12 de- 
sign cycles. At this point, the transient leading 
edge patch was removed, and the design contin- 
ued with no geometrical constraints. The 25 de- 
sign cycles required approximately 4.5 minutes 


of CPU time, while the target solution required 
approximately 2.0 minutes. This is quite effi- 
cient since the design lift coefficient converged 
to 98% of the target value in the same amount 
of time required to perform the baseline analysis 
of the target airfoil. 



0.00 0.25 0.50 0.75 1.00 
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Figure 7. Inviscid airfoil design results 
(Moo = 0.30, a = 4.00°). 


Viscous Airfoil Design 


A viscous airfoil design case was next ex- 
amined for: M^ = 0.30, Re = 2.0 ,X 10 6 , 
a = 6.00°. The target airfoil is the NAC'A 2412, 
with the flow assumed to be fully turbulent, and 
attached. The one-equation Spalart-Allmaras 
turbulence model was used. The grid dimen- 
sions for this test case were 289 X 49, with the 
y + value for the first grid point off of the surface 
being approximately 1.0. 

Figure 8 compares the results of this test 
case to the initial and target values. The SC2D 
method has again provided accurate results in 
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just 25 design cycles. The lift coefficient con- 
verged monotonically to almost 99% of the final 
value in just 12 design cycles, which would be 
quite adequate for most design problems. This 
design case required approximately 6.0 minutes 
of CPU time. This is again efficient, since the 
target solution required approximately 3.0 min- 
utes of CPU time. Except for the increased flow 
solver computational time, the viscous effects 
did not alter the stability of the SC2D design 
method. 



0.00 0.25 0.50 0.75 1.00 


x/c 

Figure 8. Viscous airfoil design results 
(M (>0 = 0.30, Re = 2.0 X 10 6 , a = 6.00°). 

High-Lift Flap Design 

The SC2D design method was initially 
developed using a two-dimensional, constant 
strength source-doublet panel method [11] as 
the flow solver. Since panel methods still play 
an important role in aerodynamic design meth- 
ods, a simple high-lift design example was for- 
mulated. Figure 9 shows the high-lift config- 
uration which was adapted from Reference 12. 
The flap is deflected 20.0° and the chord length 


is 20% of the stowed airfoil chord, ft should be 
noted that the flap gap and overlap with respect 
to the main element were arbitrarily set. 



Figure 9. Simple two-dimensional high-lift 
configuration. 

The objective of this design case was to de- 
sign the entire flap geometry, ft must be real- 
ized that an aircraft designer cannot typically 
alter the entire flap, as portions are constrained 
by the cruise wing configuration. Nevertheless, 
this test case will further demonstrate the flex- 
ibility of the new design method. The a.ngle- 
of-attaok for the test case was 0.00°, and a free 
stream Mach number of 0.20 was used in SC2D. 
The initial flap geometry was the NAC'A 0012, 
while the target flap was the RAE 5212. 



Figure 10. Incompressible high-lift flap de- 
sign results (a = 0.00°). 


The results of this test case are shown in 
Figure 10, where they are plotted in the unro- 
tated reference frame. The new design method 
has once again provided accurate results. This 
is a challenging test case not only because of the 
strong coupling between the flap and main ele- 
ment, but also because the stagnation point is 
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located considerably aft of the leading edge on 
the lower surface. This case required 200 design 
cycles to converge which is attributed to the ex- 
plicit Kutta condition used in the panel method, 
which required a lower value of a (0.25) to avoid 
kinks in the trailing edge region. As a result of 
this reduced a value, the transient leading edge 
patch was not required. The final flap geome- 
try is compared to the target geometry in Figure 
ff, where the agreement is observed to be ex- 
cellent. Even though more design cycles were 
required, this test case used only 2.0 minutes of 
CPU time. This suggests that the SC2D design 
method may be a useful tool for designing or 
improving high-lift systems. Further research is 
being conducted to explore this possibility. 



get flap geometries. 

Inviscid Transonic Airfoil Design 


The final two-dimensional test case pre- 
sented is an inviscid transonic design, with the 
RAE 2822 as the target airfoil. The test case 
was formulated to produce a shock wave free 
flow held, since the terms necessary to handle 
supersonic how [3] were not added to equation 
(9). The how conditions are M (>CJ = 0.725 and 
a = 0.00°, which generates a lift coefficient of 
0.382. The grid used for this test case was sim- 
ilar to the previous Euler case. 

Figure 12 compares the design results to the 
initial and target values. The SC2D method has 


again done an excellent job of replicating the 
target pressure distribution. The design airfoil 
geometry is identical to within plotting accu- 
racy of the target airfoil, which is quite encour- 
aging. This case however did require nearly 50 
design cycles to converge, which is attributed 
to the small sonic region present on the upper 
surface of the target airfoil. In such a region a 
slope based term [3] would be more appropriate 
than equation (9), and could easily be added. 



0.00 0.25 0.50 0.75 1.00 


x/c 

Figure 12. Inviscid transonic airfoil design 
results (M (>Cj = 0.725, a = 0.00°). 


This design case required 5.0 minutes of 
CPU time, whereas the target solution required 
approximately 1.0 minute. This is still consid- 
ered an efficient design case, particularly when 
the accuracy of the results is considered. This 
test case indicates that the SC2D method may 
be a viable tool for designing transonic airfoils, 
including advanced leading edge concepts [13]. 
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Inviscid Wing Design 

The final test case is a three-dimensional in- 
viscid wing design. The wing has an aspect- 
ratio of 7.0, taper ratio of 0.70, leading edge 
sweep angle of 25°, and no twist. The initial 
wing is comprised of the NAC'A 00f2 airfoil sec- 
tion. The target wing is generated by replacing 
the root section with the RAE 52 f 2 airfoil as 
shown in Figure 13. The flow conditions for 
the design case are: M (>CJ = 0.20 , a = 6.00°. 
The C-0 grids used to represent the wings are 
based on the grid refinement studies presented 
in reference [10]. The grid has dimensions of 
193 X 33 X 49 in the stream wise, normal, and 
span wise directions, respectively. 



Figure 13. Planform view of target wing. 

The results of the wing root section design, 
'2y/b = 0.00, are presented in Figure 14. The 
SC2D method has again performed well, con- 
verging to the target pressure distribution in 
only 25 design cycles, although the geometrical 
changes are significant. The lift coefficient again 
converged monotonically, and reached 99% of 
the final value after only 12 design cycles, ft is 
important to note that this was the only design 
station on the wing, and that existing options 
in CD1SC were used to interpolate the design 
changes between this station and the fixed outer 
portion of the wing at '2y/b = 0.50. This de- 
sign case required 40.0 CPU minutes, while the 


target solution required 15.0 minutes. This is 
again considered to be efficient since the design 
converged to 99% of the target lift coefficient in 
almost the amount of time required to perform 
the baseline analysis of the target geometry. In 
addition, the new design method behaved in a 
manner consistent with the previous low-speed 
two-dimensional design cases discussed above. 



0.00 0.25 0.50 0.75 1.00 


x/c 

Figure 14. Inviscid wing design results 
(M (>0 = 0.20, a = 6.00°). 

Conclusions 

A new streamline curvature design method 
was developed which is based on a new method 
of modeling the streamline curvature variation 
normal to the surface. The surface coordinates 
are allowed to move normal to the surface, and 
the method is valid in the leading edge region 
of airfoils. The new design method was incor- 
porated into the CDISC design method. Two- 
dimensional and three-dimensional test cases 
were presented which demonstrated the valid- 
ity of the new design method. 
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The new SC2D design method was shown 
to accurately replicate the target surface pres- 
sure distributions of the design test cases. The 
resulting design geometries were found to be 
within plotting accuracy of the target geome- 
tries, and were obtained with no geometrical 
constraints. The performance of the new design 
method, particularly in the leading edge region, 
suggests that the method may be a viable tool 
for designing both low-speed and transonic lead- 
ing edges. 

The computational efficiency of the new 
SC2D design method has been demonstrated. 
The design cases were found to converge in a 
monotonic fashion, with the lift coefficient con- 
verging to at least 98% of the target value in 
a similar amount of time required to perform 
the baseline analysis of the target configuration. 
This is quite encouraging given that most of 
the design cases required significant geometri- 
cal changes. 
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